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1 Introduction 

In a quantum system the related commutation rules determine the character of basic operators 
involved in the model, e.g. boson, fermion, spin, anyon etc. When such quantum systems are 
integrable the commutation rules are determined by the underlying Yang-Baxter (YB) algebra 
obtained from the quantum Yang-Baxter equation (QYBE) 

R(X - n)Lj(X) ® Lj(fi) = (I ® Lj(fi))(Lj(X) <8> I)R(X — fi), j = l,2,...,N (1.1) 

together with the ultralocality condition Lj(X) ® L^fi) = (I <8) Lf s (fj,))(Lj(X) ® I), at k ^ j, 
which are at the same time sufficient for the integrability of the system. Here Lj(X) is the 
representative Lax operator of the lattice (or discretized) quantum integrable systems (QIS) at 
each lattice site j = 1, 2, . . . , TV and R(X — fj,) is the quantum i?-matrix (c-number matrix) which 
determines the structure constants of the YB algebra. For rational solution of the i?-matrix 
one usually gets the spin algebra or its bosonic realization, while the trigonometric solution 
yields (/-spin and (/-bosons, which are comparatively new entries in the field of quantum physics, 
discovered mainly from the study of integrable systems [H HI [3] . 

The concept of (/-boson was formally introduced as (/-parameter deformation of a standard 
boson, or through two-mode realization of quantum algebra [HEIE]. Subsequently, the fin- 
ding of more application oriented single mode (/-bosonic realization of quantum algebra through 
(/-deformation of Holstein-Primakov transformation [9] and linking it to QIS [H)] was 
a significant step. However, it seems that this initial but important step has not been pursued 
with enough zest to explore the possibility of appearance of (/-boson in other QIS and its app- 
lications in exactly solvable physical models. It is true that there were attempts to construct 
interesting (/-bosonic and (/-spin models in various fields and to study their physical effects 
[3 El CEU [HI [I3j [TU CE]. However most of such models could not be related to the underlying 
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YB algebra of a QIS permitting exact treatment by the Bethe ansatz (BA) method. Our aim 
here therefore, is mainly to identify the (/-bosonic mode in a variety of QIS, which is deeply 
rooted in discrete or lattice regularized models allowing BA solutions. We discover also the 
appearance of such modes in some new QIS, e.g. physically relevant integrable matter-radiation 
models, a coupled double-mode (/-bosonic model leading to a two-component derivative NLS 
model etc. and furthermore by introducing a q-anyon like concept find its application in another 
new exactly solvable derivative-^ function ID anyon gas. 



2 g-boson and qr-spin 

We review in brief the basic formulation of the (/-boson and its relation to various other objects 
like (/-spin, standard boson, canonical variables etc. (/-boson can be defined through a deforma- 
tion of the bosonic commutation relation (CR): 

b q b\-q~%b q = q N , [b q ,N] = b q , [b\,N] = -b\, (2.1) 

which is supposed to be invariant under a reflection q — ► q" 1 . Therefore combining such two 
relations we can define the (/-boson in a more symmetric form 

h ^ ] = cos(a(2N + l)) ^ KN]=bq) [b l,N] = -bl (2.2) 



where q = e 2ta . At q — ► 1 the (/-boson goes into a standard boson with both the above relations 
reducing to the bosonic CR: 

[b,tf] = l, [b,N]=b, [tf,N]=-tf. (2.3) 

One can find also a direct mapping between (/-deformed and undeformed bosons as 

b q = bf(N), b\ = f(N)b\ /(JV)=^MiV N = b\ (2.4) 

where [x] q = q q Z q ~i ■ Using these relations and the usual representation of bosonic operators 
it is easy therefore to construct a (/-bosonic representation [JJ [5l [6]. Initially the (/-boson was 
introduced as a two-mode realization of the g-spin (g-deformed Schwinger transformation) 

4 = ^A 2 , s- = b\ 2 b ql , s 3 = -(N 1 -N 2 ) 
linked to quantum algebra su q {2) 

[s+,s-}=e[2s%, [ S 3 , S ±]=± S ±, (2.5) 

for e = 1. e = — 1 on the other hand gives noncompact algebra su q (l, 1). Subsequently, a single 
mode realization was found through g-Holstein-Primakov (g-HP) transformation El E] 

s s = stN, s+ =g±(N,q)b q , s~ =b\g±(N,q), 9±(N, q) = [2s =F N] q , (2.6) 

for su q (2) and su q (l, 1), respectively, with s-being a spin parameter. One can find also mapping 
from q-boson to canonical variables: [u,p] = i as 

b q = e-*/ c («), b\ = f c (uyv, f c {u) = [c- u] q (2.7) 

with c = const. 
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3 g-deformed Yang— Baxter algebra 

The Yang-Baxter algebra underlying QIS and providing guarantee for the integrability of a dis- 
crete or lattice regularized quantum system was found to be given in the (/-deformed case by 
a general form [161 El 02] 

[S+, S~] = (M + sin (2aS 3 ) + M~ cos (2aS 3 )) 

[S 3 ,S±]=±S± [M ± ,-]=0. (3.1) 

Note that this is a novel type of deformed algebra, which apart from an usual parameter q-de- 
formation, also involves additional operator-deformation by central operators M^. Operator 
deformed quantum Yang-Baxter algebra (|3.ip is a bialgebra and exhibits noncommutative Hopf 
algebra properties with the coproduct structures 

A(5+) = 5+ fg> cfg~ sS + c^q s " ® S+, A(5~) = 5" c+<?~ SJ + cfq sS ® S~, 

A(S 3 ) = /®S 3 + S 3 ®/, A(c^) = c±®c± a = 1,2, (3.2) 

where c„ , a = 1,2 are related to central operators as M ± = ^(c^Cg ± c^cj ) and one redefines 
(/ = e m . Note that unlike deforming parameter q, deforming central operators and hence 
have nontrivial coproduct. 

The quantum Lax operator associated with this general YB algebra is 

L &e iaS3 + t^e- 53 2 sin trf" \ A 

W0 " ^ 2sinaS+ £c+ e —* 3 + r 1 ^**" J ' C ( } 

This ancestor model [161 [T71 [T8] seems to be general enough to generate all quantum integrable 
models with (2 x 2) Lax operator associated with the trigonometric R tvls matrix 

R\\ = rH = s in(A + a), R{% = F%{ = sin A, R\\ = R\\ = sin a. (3.4) 



Therefore identifying the situations when the YB algebra f)3. 1 j) goes to different known algebras 
at different values of the deforming central operators, we can construct from the corresponding 
reductions of the Lax operator f|3.3|) a series of quantum integrable models solvable by the Bethe 
ansatz. 

For example we may observe that in contrast to the well known (/-deformed algebras (|2.2p . 
(|2.5|) . algebra (|3.ip has two different terms in the RHS of its main relation: the first one (sine- 
term) is similar to the quantum-spin algebra (|2,5p . while the second term ( cosine-term) to the 
(/-boson (|2,2p . Indeed both these (/-deformed algebras can be obtained as subalgebras of (|3.ip . 
Notice that M + = 1, M~ = yields quantum algebra sl q (2), while M + = — 1, M~ = gives 
the noncompact algebra sl q (l, 1) as defined in (|2.5p . The Lax operators of the related integrable 
systems can be obtained directly from the general form (|3.3p . by using a compatible choice 
cf = £2 = 1, cj~ = = e, e = ±1 as 

Z^ + e^e-^ 3 2 sin as- \ 

Lqs{Q ~ { 2sina S + ie~^ 3 + e^e™ 3 ) ' (3 ' 5) 

involving (/-spin belonging to sl q {2) and sl q (l, 1), for e = ±1, respectively. 

On the other hand for the particular values of the deforming operators: M + = sin a, M~ = 
icosa, one directly obtains the (/-boson (|2.2p as a subalgebra of (|3.ip with a compatible choice 

= = 1, = — iq, c-2 = -■ However note that such a choice of deforming operators, 
as seen from (|3.2p . seemingly does not respect the usual coproduct structure and hence the 
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corresponding algebra representing g-boson exhibits the known difficulty in its formulation of 
the Hopf algebra. We get the g-bosonic Lax operator from the reduction of (|3.3|) as 

T t^-( & laN ~ »e-l e -*«(^-l) Kb\ \ 

with a direct mapping: 

S, + = ^&„ S- = ^b\, S 3 = N, Q = q-q^ (3.7) 
and k 2 = — • Interestingly if on the other hand we choose 

= &2 = with c± = £2 = (3.8) 

or at least any one of is 0, by representing 



S+ = Kb q qf, S-=kqTbl S 3 = N, k 2 = ^-, (3.9) 



we can get a direct realization to the g-bosonic relations (|2.ip as a subalgebra of (|3.ip . We 
construct also the corresponding Lax operator from the related reduction of f|3.3|) . We see below 
that such g-bosonic models are intimately related to a series of interesting QIS including the 
quantum lattice Liouville model. 



4 g-bosons in discrete QIS 

We identify and present here mostly quantum integrable g-bosonic models which can be linked 
directly to Yang-Baxter algebra and solved exactly through algebraic Bethe ansatz. 



4.1 q-Bose gas models 

i) A quantum integrable model involving g-bosons constructed on a lattice was proposed in |lUj . 
which can be defined simply through Lax operator (|3,5p involving g-spin operators (s^,s 3 ) G 

sl q (2). g-spin operators are in turn realized in single-mode g-boson b q , b\ using g-HP transfor- 
mation (12. 6p . This lattice g-Bose gas model associated with (|3.5p and R tng exactly satisfies 
the QYBE by construction (since a particular case of ancestor model (|3.3p ) and therefore the 
related eigenvalue problem can be solved exactly through algebraic Bethe ansatz (ABA). In 
analogy with the exact lattice model of nonlinear Schrodinger (NLS) equation [19] the g-Bose 
gas model on the lattice can be shown to have localized Hamiltonian [10], as obtained from the 
lattice NLS. However at the continuum limit the g-bosonic mode turns into bosonic field and 
the model goes to the same NLS field model. 

ii) A g-harmonic oscillator model H q = ^-(\_N\ q -\-[N +\\ q ) was considered in [21], motivated 
by the noncommutative quantum field theory, in analogy with that of the harmonic oscillator. 
Though this model does not belong to quantum integrable systems, it allows to calculate explicit 
result on different thermodynamic properties like free energy, specific heat, entropy etc. 



4.2 Quantum lattice Liouville model as q-bosonic model 

Quantum integrable exact lattice version of the Liouville model may be constructed as reduction 
of the above ancestor model, when the deforming operators are chosen as (|3.8p as in the case of 
g-bosons and generators of algebra (|3.ip are realized as 



S* = iu, S+ = e-^f(u), S q =f(u)e^, 



(4.1) 
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where f(u) 2 = 1 — e 2 q~ 2lu+l . The Lax operator of this lattice Liouville model can be obtained 
directly from (|3.3|) and therefore its exact quantum integr ability is guaranteed. If one relates 
the parameter with the lattice constant: e = iA, then it is easy to show that the model goes to 
the Liouville QFT model at the continuum limit with the field Lax operator recovered from the 
present discrete Liouville Lax operator. 

On the other hand, it is interesting to observe that, the g-boson realization (|3.9p . using its 
mapping (|2.7p through the canonical variables turns unexpectedly to realization (|4.ip linked 
to the lattice Liouville model. Therefore we may conclude that the lattice Liouville model is 
intimately related to and can be represented by the (/-bosonic model. 



4.3 Quantum derivative NLS through q-boson model 



Reduction of the ancestor model to Lax operator (13. 6ft expressed directly in (/-bosons gives 
another quantum integrable deformed Bose gas model on the lattice. L-operator (|3.6|) together 
with R g -matrix (|3.4p satisfy the QYBE by construction and allow exact solution of this g-bo- 
sonic model by algebraic AB (ABA). Using mapping (|2.4f) the g-boson defined on a lattice can 
be linked to a bosonic operator with commutation relation [^j,^)t] = -^i in the form 

It is interesting to note that at the continuum limit A — > 0, when the lattice boson tpi produces 
the bosonic field tpi — > y/Aip(x), the (/-boson through realization (|3.7p with the choice 

c i = 4 = 1, c i = - J*?. c 2 = j- ( 4 - 3 ) 

goes to a quantum integrable derivative NLS field model [20], given by the equation iipf = 
ipxx — i>)^x- The lattice Lax operator (|3.6|) with mapping (|4.2|) . (|4.3p reduces at A — > to 
L = 1 + AC/ DNLS , where 



U 



DNLS 




is the Lax operator of the DNLS field model [20J. Note that the quantum DNLS Hamiltonian 
in the iV-particle sector is equivalent to the interacting Bose gas with derivative <5-function 
potential [45j. We introduce in Section 7 a (/-deformed anyon on the lattice and subsequently 
construct an anyonic gas model interacting through <5'-function, which is solvable by exact Bethe 
ansatz. 



4.4 Coupled q-bosons and multi-component quantum integrable 
derivative NLS model 

Note that the Yang-Baxter algebra (|3.ip is invariant under the exchange of «-> c\f, since 
do not change under such transformation. However, this is not the case with the associated 
Lax operator (|3.3p . which is transformed as L(£) — ► L _1 (|). Therefore, using such a symmetry 
of algebra (|3.ip we can introduce another Lax operator from (|3.6|) with a similar but different 
(/-boson. Fusing these two Lax operators we can construct a novel quantum integrable two-mode 
coupled g-bosonic lattice model given by the Lax operator L 2 i b {£) = Li b (£,b ql )(Li b (±,bq 2 )y 1 
with explicit expression of its matrix elements as [22] 



6 



A. Kundu 



Lit = n Q (q-»% - 4^ 3+1 ) +f0{l*-" a ) " 4 6 ^ 1+1 ) ' (44) 

with = (-L^i )^ and L^ 6 = (-^12 & )^- By combining N + number of g-boson of the first kind 
and N~ number from the second kind one can construct also a multi-mode generalization of 
this g-bosonic model. The mutually commuting conserved operators Ck, k = 0, ±1, ±2, .... in 
the simplest case of two-mode lattice g-bosonic model may be generated by 

through expansion in spectral parameter: r(£) = ^Cfc£ fc . This long-range interacting lattice 

h 

two-mode g-bosonic model is a new quantum integrable model and should be solved through 
algebraic Bethe ansatz for exact result. The possibility of obtaining its local Hamiltonian with 
few nearest neighbor interactions by applying the method developed in [19] and used for g-Bose 
gas in [10] is also an interesting problem to explore. Considering A to be the lattice constant, 
we can go to the continuum limit by taking A — > 0, where the two-mode g-bosonic lattice Lax 
operator (fOj) turns into L 2qb (£) ->• I + iAU M ™(£,). Interestingly, Lax operator £/ M ™(£) thus 
obtained is the same field Lax operator associated with the bosonic massive Thirring model as 
well as with a two-component derivative NLS model for different choices for the Hamiltonian. 
The elements of U M ™ {£) are given bv [231120] 

jjMTM = _jjWVM = 1 Q. _ ^ + _ 

We can evaluate the commuting conserved quantities of this integrable field model at the classical 
limit (putting k± — > 1) using the Riccati equation: 

U X2 (j^j + 2U n u + u 2 = U 12 U 2l . 

We solve this equation by expanding u = ^ C 2 kC 2k i n spectral parameter at £ — ► and at 

k 

£ — ► oo, obtaining the conserved charges Co, C2, C4, . . . and C_o, C_2, C_4, . . ., respectively. In 
explicit form they are given as 

Co = |0 2 | 2 , 

C 2 = ^ x -^ 2 \ 2 \^\ 2 -2^ 2x , 

Ca = _c 2x - 4</> 2 <^i| 2 + ^ ~ + m 2 + l^ 2 | 2 )(l - 2C 2 ) + <f>f^ 

and 

C- O = -|0i| 2 , 

C_ 2 = 0^ 2 -4|0 2 | 2 |0i| 2 + 2^^, 

c_ 4 = c_ 2a; + H^M 2 - ^ (Vl + ^2) - (I0il 2 + l^ 2 | 2 )(i + 2C_ 2 ) - ^ 2 |^. 

It is interesting to observe that two different combinations of the above conserved operators 
of this integrable system can generate two important field models, which can also be raised 
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to the quantum level as quantum integrable systems solvable by the algebraic Bethe ansatz. 
For example the field Hamiltonian given by -Hmtm = J dx(C2 + C_2) and the momentum 
P = J dx{p2 — C-2), with canonical PB relations {4>k(x), <j)\ (y)}pB = i5ki5(x — y) would yield 
the relativistic massive Thirring model (MTM), rather its bosonic version. On the other hand 
the higher Hamiltonian i?2DNLS = / dx(C4 — C-4) would generate an integrable 2-component 
DNLS field model. At quantum level it would be an integrable QFT model, an exact lattice 
regularized version of which we have constructed already through two-mode (/-bosonic model. 
It is desirable to investigate this new quantum integrable model both at the discrete and the 
continuum limit, through exact Bethe ansatz formalism. 

4.5 q-boson and Ablowitz— Ladik model 

The celebrated Ablowitz-Ladik (AL) model was discovered as a discretized version of the NLS 
field model [23] much before the discovery of (/-boson. However surprisingly one identifies 
that the underlying YB algebra of the quantum integrable AL model is actually given by q- 
bosons. We observe first that twisting transformation R^j(X) — > e id ^~ k ^R^j(X), can give a new 
i?-matrix solution, where 9 is some free parameter. Using this twisted i? trig -matrix one gets 
a ^-deformation of the quantum algebra (|3.ip . where the commutator is deformed as pSgS~ — 
p~ 1 S q 'S^, p = e l9 , (similar to p, (/-deformation) [THJ [TTJ [18] . As a result one gets a ^-deformed 
(/-boson. For a special value 9 = —a and with a redefinition of (/-boson as b q = 6 <? g i_ ^"0 1 / 2 
and bq = q~^bq$l l l 2 , one gets another form of (/-boson with algebra q 2 b q b\ — b\b q = q 2 — 1 
as introduced by Macfarlane 0]. The Lax operator of the corresponding integrable model as 
reduced from (|3.3[) is given simply by 



with h = 1 — q~ 2 , which turns out to be the lattice model proposed by Ablowitz and Ladik [23] 
and its quantum generalization discussed in [25]. 

We present another novel class of important models involving (/-bosons, e.g. (/-deformed 
matter-radiation models, in the next separate section. 

5 Integrable matter-radiation models with qr-bosons 

Matter-radiation (MR) models represented by atoms interacting with radiation are usually de- 
scribed by two-level spin operators interacting with single mode boson. The well known and 
simplest models of this type are the Jaynes-Cummings (JC) [26], and the Buck-Sukumar (BS) 
model [27] . The basic physics underlying a variety of important phenomena in interacting MR 
systems, like those in quantum optics induced by resonance interaction between atom and a 
quantized laser field, in cavity QED, in trapped ion interacting with its center of mass motion 
irradiated by a laser beam [28] etc., seems to be nicely captured by such simple models. Many 
theoretical predictions based on these models, like vacuum Rabi splitting, Rabi oscillation and 
its quantum collapse and revival etc. have been verified in maser and laser experiments. How- 
ever, for describing physical situations more accurately generalizations of these basic models, 
like (/-deformed BS and JC model [3 El [29J [30] , (/-boson model interacting with (/-spins [31] or its 
classical variant [32], trapped ion (TI) with nonlinear coupling, multi-atom models [331 [M] etc. 
have been proposed. However most of the above generalizations, except a few [31], go beyond 
exact solutions and quantum integrability, especially for interacting multi-atom models. Various 
(/-bosonic models and quantum group symmetries were discussed in [351 l36l [37] . 
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We propose new quantum integrable generalization of MR models involving g-bosons and 
(/-spins. A general approach for constructing integrable MR model has been reported in our 
earlier work [38]. The physical significance of g-deformed MR models in comparison with their 
undeformed counterparts is given by their stronger nonlinear interactions between atomic ex- 
citations and the radiation mode, as well as the possible presence of nonlocal and asymmetric 
interactions between atoms due to nontrivial coproduct structure of (/-spins. For constructing 
these models in a unified way and for ensuring their quantum integrability we may start from 
the integrable ancestor model (13. 3h or more precisely by forming their Lax operators through 
a combination of (|3.5p involving (/-spins and (|3.6[) linked to the (/-boson. The general form of 
such integrable (/-deformed MR models may be given by the Lax operator 

L = £ a nc(£> S q , S 3 )L qs (£, s q , s 3 ), (5-1) 

where -L anc (£, S^, S 3 ) is given by (|3.3|) and L gs (^,s^,s 3 ) by (|3.5p with e = +1 . Therefore by 
construction it would exactly satisfy the QYBE with tr (L) = ^ C n £ n defining the mutually 

n 

commuting conserved operators C n , n = 0,1,2, ... of this quantum integrable system. The 
Hamiltonian of the model given by Cq takes the form: 

H qM R = H d + (s+S~ + s-S+)r), H d = c + cosa(S 3 -s 3 ) + ic~ sin a(S 3 -s 3 ), (5.2) 

where rj = 4 sin 2 a, and c are related to the values of the deforming operators appearing in 
the Lax operators. This is a quantum integrable system with another commuting conserved 
operator C2 = S 3 + s 3 . Note that the quantum-spin operator can be expressed through iV-num- 
ber of spin-^ operators using the coproduct structure defined in the tensor product of N vector 
spaces as 

s - E 4 E of 3 

This describes an asymmetric as well as nonlocal interaction between the two-level atoms me- 
diated by the (/-bosonic radiation mode. 

Different reductions of this generalized g-deformed MR model performed through possible 
realizations of the generators S^, S 3 of the ancestor algebra, as we have listed above, can 
generate different physically relevant MR models in a unified way, as we describe below. 

5.1 Integrable q-deformed Jaynes-Cummings model 

This model is constructed as a quantum integrable system of interacting (/-spins with (/-boson 
representing atoms interacting nonlinear ly with radiation as well as among themselves. The Lax 
operator of this integrable model can be constructed as a reduction of (|5.ip by replacing the 
ancestor model by (13. 6p using the (/-bosonic realization (13. 7p . Consequently the Hamiltonian of 
the model can be constructed from (|5.2p as 

H q jc = c sin a(N - s 3 + to) + {b\s+ + s-b q )r], (5.4) 

where c, rj, u are constant parameters (dependent on a), adjusted to simplify the expression. 
Notice that at q — ► 1 limit, when — > s^, b q — > b, the g-deformed Jaynes-Cummings (gJC) 
model (|5.4p goes to the integrable multi-atom Jaynes-Cummings model (at a 2 order). In general 
this model describes nonlinear and nonlocal interactions between atoms and the radiation mode, 
since can be expressed through iV-number of two-level atoms in a nonlocal way as (15. 3D and 

the radiation represented by the g-bosonic mode b q , b\ is expressed through standard boson b, b^ 
in a nonlinear way through the mapping (|2.4p . 
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Another similar multi-atom qJC model with explicit interatomic interactions may be given 
by defining the Lax operator as 

N 

L = L qh J[L<il, 

3 

(j) 1 

where Vixz is the Lax operator of anisotropic xxz spin- ^ chain, which represents here the N 
interacting two-level atoms. The corresponding Hamiltonian of this model would include more 
complicated matter-radiation interactions as well as explicit atom-atom interactions. 

Both the above qJC models are novel quantum integrable models solvable exactly by algebraic 
Bethe ansatz. Detailed analysis of these models with possible physical importance needs further 
pursuing. 

5.2 Integrable q-deformed Buck Sukumar model 

This quantum integrable model also describes matter-radiation interaction through (/-spins and 
a (/-boson interacting in a stronger nonlinear way. The idea of construction is to start from 
the Lax operator (|5.ip and take the realization of the ancestor model through generators of 
the noncompact quantum algebra su q (l, 1) as given by (|3.5p for e = —1. Subsequently such 
generators are realized through (/-boson via (/-Holstein-Primakov transformation: 

to represent the radiation mode. Consequently the Hamiltonian of this (/-deformed Buck- 
Sukumar ((/BS) model which can be reduced from the general (/MR model (15. 2j) takes the form 

H qBS = csma(N - s 3 + u) + rj (j[N] q bls q + s+b q J\N]?) . 

The g-deformed BS model gives an integrable version of an earlier model [8], when (/-spin 
operator is replaced by <r-matrices by taking spin-i representation. 

This novel quantum integrable (/BS model can be exactly solved using Bethe ansatz, which 
should produce physically interesting result generalizing that of the well known BS model, and 
therefore it deserves detailed investigation. At q — > 0, as is clearly seen, the (/BS model goes to 
the integrable multi-atom BS model. 

Note that the quantum integrable model of (/-bosons interacting with (/-spins proposed in [31] 
is similar in spirit to the present model, where however (/-bosons were introduced not directly 
but as realization of su q {2) through (/-HP transformation. 

5.3 Integrable (/-deformed trapped-ion model 

Though this model does not involve (/-bosons directly, we present it here since this novel quantum 
integrable MR model belongs to the same trigonometric class associated with the quantum 
R g -matrix and can be obtained again from the (/-deformed MR model (|5.2p under suitable 
realization. 

Quantum Yang-Baxter algebra (|3,ip under reduction = 0, when both = 0, simp- 
lifies to 

[S q , S q ] = 0, [S 3 , S q ] = zLSg 1 , 
allowing realization through canonical operators: [x,p] = i, as 
Sf = e^ iX , S 3 = p. 
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Note that the related reduction of (|3.3j) yields the Lax operator of the quantum integrable 
relativistic Toda chain [411 1421 1^1- Interestingly the same realization can reproduce from the 
general MR Hamiltonian (|5.2p a (/-deformation of the trapped ion ((/TI) model given by 

H q Ti = c+ cos a(p — s 3 + uj) + c- sin a(p — s 3 + uj) + r]{e~ lx s~ + e lx st) 

with highly nonlinear coupling between the atomic excitation and the vibration of the center 
of mass motion of the trapped ion, described by displacement x. The system also has another 
conserved operator C\ = p + s 3 . Usually for achieving exact solution such nonlinear oscil- 
lations are linearized through several approximations like Dicke approximation, rotating wave 
approximation etc. [28] . The present model however is solvable in principle with full exponential 
nonlinearity and without any approximation. 



6 Exact solution through algebraic Bethe ansatz 

Almost all (/-deformed models presented here similarly to their unified construction allow exact 
ABA solutions also in a unified and vastly in a model-independent way. 

Note that from the local QYBE taking the tensor product T(A) = Y\ L>j(\) one can go to its 

j 

global form: 

R{\- n)T(\)®T(p) = {I®T(fi))(T(\)®I)R(\- fi), j = l,2,...,N (6.1) 

reflecting the Hopf algebra structure of the underlying YB algebra. Taking trace from both 
the sides of ()6.1j) and defining r(A) = trT(A) we get [r(A),r(^)] = and expanding further in 
spectral parameter: r(A) = ^ C n X n derive finally the quantum integrability condition for the 

n 

conserved operators: [C n ,C m ] = 0. 

ABA formalism aims to solve exactly the eigenvalue problem for all conserved operators si- 
multaneously. The diagonal entries r(A) = Tn(X) + T22(A) produce the conserved operators, 
while the off-diagonal elements T2i(A) = -B(A) and Ti2(A) = C(A) act like creation and annihila- 
tion operators of pseudoparticles. The M-particle state is defined as \M) b = B(X\) ■ ■ ■ B(Xm)\0) 
and the pseudovacuum |0) is defined through C(A)|0) = 0. The basic idea of algebraic BA [43] 
is to find the eigenvalue solution: t(X)\M)b = A(A, {A a })|M)#, for which diagonal elements 
Tjj(A), i = 1, 2 are pushed through the string of -B(A a )'s toward |0), using the commutation 
relations obtainable from the QYBE (|6.ip . Considering further the actions Tn(A)|0) = a(A)|0), 
722(A)|0) = /3(A) 1 0) , one arrives finally at the eigenvalue expression 

M M 

A(A, {A a }) = a(A) J] /(A - A a ) + /3(A) J] /(A a - A), (6.2) 

a=l a=l 

where /(A) is defined through the elements of the i? t rig-matrix as &in s ^"^ . Expanding A(A, {A a }) 
in powers of A we obtain the eigenvalues for all conserved operators including the Hamiltonian, 
where the rapidity parameters {A a } involved can be determined from the Bethe equations 

W7)=^ a J{X^-y o = 1, 2, . . . , Af, (6.3) 

which follow in turn from the requirement of \M)b to be an eigenvector. Returning to our 
models we find that, the major parts in key algebraic BA relations (|6,2p and (|6.3p . described 
by i?-matrix elements /(A) is model-independent and hence same for all g-bosonic and (/-spin 
models belonging to the trigonometric class. 
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The only model-dependent parts in these equations, expressed through a(A) and /3(A) are 
determined from the Lax operator and therefore can be constructed in a unified way starting 
from (|3.3p as 

«(A) =4tq m +cl-q- m , /3(A) =c+Cq- m + cdq m , (6.4) 

where m = — (0|S 3 |0) denotes the (/-spin projection. 

Since the q-bosonic or (/-spin models are obtained as various reductions of this ancestor model 
one can obtain the corresponding exact result using the proper reduction of (|6.4|) . Note however 
that for some models like relativistic Toda chain, (/-deformed trapped ion model etc. since there 
is no easy pseudovacuum construction, the standard ABA method is not applicable to them and 
more generalized functional BA has to be implemented [41) . 



7 qr-deformed anyon and 8' anyon gas 

We have found in Section 4.3, that a quantum integrable (/-bosonic model on a discrete lat- 
tice goes to an integrable derivative NLS quantum field model involving bosonic field opera- 
tors ip(x), ift(x). The Hamiltonian of this integrable QFT model may be given by 

H = [ dxtylfa + %ntfil>(^ x - , (7.1) 



which at the A r -particle sector can be shown to be equivalent to a derivative (5-function Bose gas 
model 

N 

H N = -J2 d 2 Xk +iK^ S(x k - xi) (d Xk + d Xl ) (7.2) 

k (k,l) 



which is exactly solvable by the coordinate Bethe ansatz 

Defining a new notion of g-anyon on the lattice we propose a similar model of Id anyon gas 
interacting through derivative 5-function and show that the model is also exactly solvable by 
coordinate BA. Let us consider a lattice of iV-sites and let (/-bosons: &„,-, j = 1,2,... ,N 

satisfying commutation relations (|2.2p for the same site: [b q j, b qk ] = <5jfc — — ■> etc., while all 
operators commute at different sites with j ^ k. We define another set of nonlocal operators as 

W E N h . , -id £ N k 

A qj = e k=i b qj , ^ qj = b\ j e (7.3) 

and easily check that at the same site these new operators behave exactly like (/-bosons with 
relations 



/it i cosa(2Nj + l) 
[Agj, A A = , [A q j, A q j\ = 0, etc. 



*v>-gp- cosa > W^J-", — (7-4) 

However for different sites j > I, separated by any distance we get an additional phase: 

AqjA q i = e l6 AqiAgj, A q -A ql = e %B A ql A q -, A q jA^ ql = e l9 A ql A q j, etc. (7.5) 

9 = recovers obviously the q-bosonic case, while 9 = n gives anti-commutator at different sites 
similar to fermions. Coupling now g-anyon to (/-boson mapping (|7.3p with that from (/-boson to 
boson (|2.4p and defining the boson- anyon transformation as 

j j 
A 3 =e * bjt A ) = b)e * , A, lh, A)a) (7.6) 
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we can derive a similar mapping from (/-anyon to anyon as 

A qj = AjfiNj), A\ 3 = f(N,)A], f(N) = j^)^ , N = A^A. (7.7) 

The commutation relations for the anyons Aa A\ at the same site are same as that of bo- 
sons ()2.3p . while at different sites they are same as (|7.5j) . 

At the continuum limit A — > 0, introducing a scaling through \J~K the (/-boson reduces to 
a bosonic field ip(x): 

b qi -> yfKip{x), Nj -» A^ f (x)^(x). 
Consequently the (/-anyon on the lattice is reduced to Id anyon field ip(x): 

A qi -» VA^(x), ATj -> A^t(a;)^(j;), 
where the anyonic field satisfies the commutation relations: 

^( Xl )^(x 2 ) = e ie < x ^^{x 2 )ij\x 1 ), 

${x x )${x 2 ) = e~ i6 < xi - x ^^(x 2 )^{xx) + 8(xi - x 2 ), etc., (7.8) 

where 

f ±1 for x > y, x < y, 
e(x-y) = \ „ (7.9) 

for x = y. 

From the introduction of the above discrete and Id anyonic and (/-anyonic operators it is clear 
that by replacing (/-bosons in (|3.6p by g-anyons one can construct formally a g-anyonic model. 
However unfortunately due to the nonlocal commutation rules (|7.5p the model turns into a nonul- 
tralocal model not satisfying the QYBE and hence becomes nonintegrable. In a similar way if 
we consider the fields in the DNLS model (17. lj) to be anyonic instead of bosonic, its quantum 
integrability would be lost immediately. 

However we find interestingly that, if instead of Bose gas we consider (]7.2[) as the anyonic gas 
interacting through J'-function the model remains exactly solvable through coordinate Bethe 
ansatz, though unlike the Bose gas the anyonic model is not a quantum integrable system 
with mutually commuting higher conserved operators. Without giving the details of this novel 
exactly solvable interacting anyon gas model, which we reserve for a separate publication, we 
mention only that the exact BA result of this model closely follow that of the Bose gas |45] with 
a redefinition of the coupling constant k involving the anyon parameter 6, as happens also in 
the case of 5- function anyon gas [39l HO] ■ 



8 Concluding remarks 

We have identified the appearance of g-bosons in quantum integrable systems, exploring from 
the well known to new models in a unified way. Some models, like Ablowitz-Ladik model, 
though well known, the underlying g-bosonic connection of this model was not obvious and 
even unexpected. In some cases involvement of (/-bosons is more direct, like in the exact lattice 
version of the quantum DNLS or the massive Thirring model and in some others, like in the q- 
Bose gas model of [10] . the connection to (/-boson is only through realization of (/-spin operators. 
Many models presented here, e.g. two-mode (/-boson model and related two-component DNLS 
quantum field model, (/-deformed JC, BS and TI models, as well as <5'-function anyon gas, are 
new quantum integrable systems with rich possibilities and deserve detailed investigation. Since 
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our objective here is to focus on various integrable models and identify the role of g-bosons in 
them, we could not concentrate on any individual model in detail, which we plan to do elsewhere. 

Finally since this article is dedicated to the memory of Vadim Kuznetsov, I would like to 
mention that Vadim 's interest was closely linked to the present investigation, related to quantum 
algebra and quantum integrable systems. Though unfortunately I did not have any personal 
interaction with Vadim, I remember having intense discussion with him through email regarding 
our common interest in formulating quantum relativistic Toda chain model [411 14"2| |4"4"| , which is 
intimately related to the q-deformed trapped ion model presented here. 
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